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Vortex Lattice Melting into Disentangled Liquid Followed by the 3D-2D Decoupling 

Transition in YBa 2 Cu 4 Os Single Crystals 

X.G. QiuO, V.V. Moshchalkov, and Y. Bruynseraede 
Laboratorium voor Vaste-Stoffysica en Magnetisme, Katholieke Umversiteit Leuven, 
Celestijnenlaan 200D, B-3001 Leuven, Belgium 

J. Karpinski 

Laboratorium fur Festkorperphysik, ETH-Hdnggerberg, CH-093, Zurich, Switzerland 

A sharp resistance drop associated with vortex lattice melting was observed in high quality 
YBa2Cu4 08 single crystals. The melting line is well described well by the anisotropic GL theory. 
Two thermally activated flux flow regions, which were separated by a crossover line B cr (T)=1406.5(l- 
T/T c )/T (T c =79.0 K, B cr in T), were observed in the vortex liquid phase. Activation energy for 
each region was obtained and the corresponding dissipation mechanism was discussed. Our results 
suggest that the vortex lattice in YBa2Cu40g single crystal melts into disentangled liquid, which 
then undergoes a 3D-2D decoupling transition. 
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The vortex dynamics in the mixed state of high 
T c superconductors (HTSCs) remains a subject of in- 
tense research because of its fundamental importance 
for physics of the vortex matter Jl],|],f|. ^ was P re ~ 
dieted theoretically [ ip|Jrj| and confirmed experimentally 
^IHJl^JllI,!!,!!,!^^ that, the vortex system under- 
goes a second order transition from a high temperature 
vortex liquid to a low temperature vortex glass or Bose 
glass in superconductors with strong disorders. In a clean 
superconductor a first order transition from vortex liquid 
to a well ordered vortex lattice occurs. 

Vortex melting in high quality high Tc superconduct- 
ing single crystals has been detected by various tech- 
niques including torsion oscillator Q, transport P,^0[, 
magnetization pd|JT^ | , neutron scattering |l3|], \x spin ro- 
tation and calorimetric measurements |L5|. Central 
to the question of flux melting is: does the vortex lattice 
melt via a first order transition? How does disorder influ- 
ence the melting transitio n frill ? It was found by Safar et 
al. U and by Kwok et al. [|l0| that the resistive transition 
of high quality YBa2Cu30*7 (Y123) single crystal showed 
a sharp discontinuity and that the transition exhibited 
a hysteric behavior. They attributed this phenomenon 
to flux melting. Although their explanation was ques- 
tioned by Jiang et al. JIt]] who found that the hysteresis 
was not necessarily an indication for a melting transi- 
tion, it could be resulted from a nonequilibrium behavior 
seen in the resistive transition. Subsequent magnetiza- 
tion |ll],[l2| and calorimetric |jl5| measurements convinc- 
ingly proved that the melting transition was indeed a first 
order one, i.e., discontinuity in the internal energy at the 
melting temperature was observed. By doing simulta- 
neous transport and magnetization measurements, Welp 
et al. |L2j and Fuchs et al. found that the temper- 
atures, where the resistance jumped, coincided with the 
melting temperature obtained from magnetization mea- 
surements on Y123 and Bi2Si2CaCu20s (Bi2212) single 



crystals, respectively. The role of disorder in the melt- 
ing transition has been checked by Kwok et al. jl0| and 
by Fendrich et al. fis|j . They found that strong pinning 
from twin boundaries or artificially introduced defects 
drove the first order melting transition into a second or- 
der glass transition. 

Another important question related to the flux melt- 
ing is how the melting propagates. Will the vortices lose 
their coherence along the c direction during the melting 
transition or will they keep the c-axis integrity and then 
lose it afterwards, at higher temperature, when thermal 
fluctuations destroy long range correlations along the c- 
direction pp|? This question has been recently pursued 
by employing the flux transformer configuration. The re- 
sults remain still controversial. Lopez et al. and Fuchs 
et al. found a coincidence of the melting transition and 
decoupling one in Y123 and Bi2212 single crystals, re- 
spectively |2l]]. However, Wan et al. and Keener et al. 
did similar measurements on Bi2212 single crystals and 
concluded that the melting took place in a two stage 
fashion §§. 

Previous studies of the vortex melting transition have 
been concentrated on Y123 and Bi2212 single crystals. 
Up to now, no detailed measurements on YBa2Cu40s 
single crystals have been reported. Meanwhile, the dissi- 
pation in the liquid state after the melting transition still 
remains to be clarified. The main purpose of this work is 
to show that in YBa 2 Cu40s single crystals, vortex lattice 
melts into liquid via a first order transition. It was also 
found that the liquid state after the melting transition 
could be separated into two regions which were governed 
by different dissipation mechanisms. Our results support 
the picture that the vortices first lose their transverse 
coherence at the melting temperature, this lose being fol- 
lowed by a disappearance of correlation in longitudinal 
direction at higher temperature. 

YBa2Cu408 (Y124) single crystals were grown by trav- 
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elling solvent floating zone method. The details for 
the single crystal growth were published elsewhere p3[ . 
These single crystals are needle-like with typical dimen- 
sions of 1.2x0.3x0.05 mm 3 . T c of the single crystals was 
about 76-78 K. Three single crystals were used for sim- 
ilar measurements. Each crystal was carefully cleaved 
to obtain optically flat surfaces with the c-axis normal 
to the sample surface. Gold wires were attached to the 
surface of crystal by using Pt epoxy. Then the crystal 
was heated in air at 100 °C for 1 hour, yielding a typical 
contact resistance below 0.5 Q. The resistance was mea- 
sured by four-probe low frequency ac lock-in technique 
with an excitation current of 0.3 mA at 17 Hz in the 
ab plane. The magnetic field was generated by a 15 T 
Oxford superconducting magnet. The direction between 
the crystalline c-axis and the magnetic field was adjusted 
by rotating the sample holder with an accuracy of 0.2 °. 
The crystal we report on here had a T c of 77.6 K and 
transition width of about 1.2 K at zero field. 




1000/T(K"') 

FIG. 1. Arrenhius plot of resistance vs. 1000/T at different 
applied magnetic fields. Solid lines are experimental data. 
From left to right: H= 0.2, 0.5, 1, 2, 3, 4, 5, 6, 7, 8, 10, 12 
T. The dashed lines are the fits with TAFF theory, h and I 
are used to identify the two TAFF regions. T cr shows how 
the decoupling transition temperature is determined. Inset: 
determination of T m from the dR/dT vs. T plot derived from 
the R-T curve for H=0.5 T. 

Fig.l shows the Arrenhius plot for the temperature de- 
pendence of the resistance measured at different fields 
applied perpendicular to the ab plane up to 12 T. A 
sharp jump in the resistance R with a magnitude of 
R/R n ~10% is clearly visible at each magnetic field. 
The magnitude of the resistive drop lies between those 
of Y123 (-20%) and Bi2212 (-0.5%), suggesting a pos- 
sible dependence of the R(T) jump on the anisotropy of 
the materials since Y124 has an anisotropy between those 



of Y123 and Bi2212. We also observed small hysteresis in 
the resistive jump upon warming up and cooling down. 
We measured resistive transition at 4 T with different 
excitation currents of 0.1 mA, 0.3 mA and 1 mA, re- 
spectively. The R-T curves measured at 0.1 mA and 0.3 
mA were essentially identical. However, that of 1 mA 
deviated from those of 0.1 mA and 0.3 mA. Similar hys- 
teric resistance jumps were earlier observed on Y123 and 
Bi2212 single crystals and were attributed to a first order 
flux melting in clean superconductors P,|l(| . To quantify 
the melting transition, we determine the melting temper- 
ature T m as the temperature where dR/dT vs. T plot 
shows a sharp peak. The inset in Fig. 1 demonstrates 
how the melting points T m are determined. The obtained 
melting temperatures are plotted against the correspond- 
ing magnetic fields in Fig. 2. The melting line can be de- 
scribed by the empirical formula B m (T)— (1-T/T c ) a with 
the exponent a between 1 and 2 |||||l(j . We find that the 
melting line is best fitted by B m (T)=31.51(l-T/T c ) 1 - 47 , 
T c =77.6 K being the transition temperature at zero field. 
This result agrees well with that obtained by Kwok et al. 
@ who found B m (T)=103(l-T/92.33) 141 on Y123 sin- 
gle crystal when the field was applied parallel to c-axis. 
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FIG. 2. The H-T phase diagram for vortices in 
YBa2Cu4 0g single crystal. Solid line and dashed line 
are best fittings to the experimentally obtained melt- 
ing line (filled circles) and 3D-2D decoupling line (open 
circles) with B,„(T)=31.5(l-T/77.6) 147 (solid line) and 
B cr (T) = 1406.5(l-T/79.0)/T (dashed line), respectively. 

To check further an interpretation of the B m (T) line 
in the framework of the vortex melting scenario, we mea- 
sured the resistive transition with the magnetic fields ap- 
plied at different angles to the ab plane, while fixing the 
field magnitude at 4 T. The angular dependence of the 
melting transition is shown in Fig. 3 where clear resis- 
tance jumps can be seen at all the angles. It is inter- 
esting to note that at small angles (6*<20°, 6 being the 
angle between ab plane and applied field), the resistive 
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transition shows an anomalous behavior above the kink. 
This anomaly was more clearly visible in our flux flow 
measurements (R-B curves) . Such kind of anomalous be- 
havior could be due to the vortex tilting instability at 
large angles and under a magnetic field of a few Tesla 
which induces a competition between intrinsic pinning 
and vortex interaction. Following the same procedure, 
we have obtained the angular dependence of the melting 
transition temperatures at 4 T as shown in the inset of 
Fig. 3. 

For a 3D vortex lattice, the melting transition occurs 
when the shear modulus c^e goes to zero, which is de- 
termined by using a phenomenological Lindermann cri- 
terium, i.e., when the mean square root amplitude \J (it 2 ) 
of the displacement of the vortex from its equilibrium 
position is larger than a certain portion of the lattice 
constant, -J (u 2 ) > c^ao where cl is the Lindermann 
number, eto being the vortex spacing. According to the 
theory for an anisotropic superconductor, T m (9) has such 
a relationship as || 
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FIG. 3. Angular dependence of the resistive transition at 
a fixed magnetic field of 4 T. From left to right: 6=90, 70, 
60, 50, 40, 30, 20, 10, 5, 2 °. Inset: Angular dependence 
of the melting transition in YBa2Cu4 0s, solid line is fitting 
by using anisotropic GL theory Eq. (1) with a Lindermann 
number cl=0.14. 



fc B T„ 



47T 2 A 2 ;6 (T ro )i?V2 £ l/4(0) ' 



(1) 



where e{6) = cos 2 (6>) +7 2 sin 2 ((9), $ = 2.07 x 10- 7 G 
cm 2 is the flux quantum, A a fc is the penetration depth in 
the ab plane and 7 = A c / X a b the anisotropy parameter. 
A best fit to the obtained data by using Eq. (0) with 
X ab = 2000 A (Ref. 24) is achieved when 7 =13.6, c L = 
0.14, which is shown as the solid line in the inset of Fig. 3. 
Clearly a quite satisfactory result has thus been obtained. 
The cl value agrees very well with those reported by 
Safar et al. || and by Kwok et al. M|. Considering 
the crystalline structures of Y123 (7=7.7) p), Y124 and 



Bi2212 (7=50-170) Q, we think that the anisotropy of 
7=13.6 is a reasonable value. 

As mentioned in the introductory part, in a clean lay- 
ered superconductor, the melting line separates the vor- 
tex lattice phase at low temperatures from the vortex 
liquid phase above B m (T). Depending on the tempera- 
ture, magnetic field and strength of disorder, the vortex 
can be melted into a disentangled 3D line liquid or an 
entangled liquid flTp6|]. It is important to note that the 
melting transition is not necessarily a depinning transi- 
tion, which means that vortex liquid can be pinned. The 
pinning force can arise from the residual pinning or from 
viscosity of the vortices. Due to the different dimension- 
ality and vortex configurations involved in these various 
liquid states, the dissipation mechanisms should also be 
quite different. A useful information about the liquid 
state can be extracted by carefully analyzing the activa- 
tion behavior of vortex liquid. 
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FIG. 4. Activation energies derived from the slopes in 
Fig.l. Filled circle and open circles are experimental data. 
Solid line and dashed line are the corresponding best fits. 

From Fig. 1 it is clearly seen that above the melt- 
ing temperature, i.e. in the liquid state, there are 
two distinct parts where InR shows a linear dependence 
upon 1/T with different slopes. This kind of behav- 



ior (p 



Poe 



-U/kT 



) is typical for thermally activated 



flux flow (TAFF) with an activation energy U(B,T) = 
U(B,0)(l-T/T c ) Such two separate TAFF regions 
with different activation energies have previously been 
observed by us on oxygen deficient Y123 thin films p7| . 

Using the same notations as in Ref. 27, we define the 
TAFF regions as region I and region h as shown in Fig.l. 
For each TAFF region, we extract the activation energies 
from the slopes of the Arrenhius plot. The obtained re- 
sults are shown in Fig. 4. We find that the activation 
energy can be best fitted by U l = -84.3 + 497.5/ \[B and 
Uq = 550— 423.6 lni?, respectively. Moreover, from a low 
temperature extrapolation of TAFF fit to region h and 
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a high temperature one for region I, we obtain the tem- 
peratures T cr (B) which define the crossover from region 
I to h. 

It was argued []T|,p6| that the Abrikosov lattice could 
melt into a disentangled liquid under certain circum- 
stances. As discussed by Geshkenbein et al. and by 
Vinokur et al., the dissipation in a disentangled liquid 
can be developed via plastic deformations of the vortices 
p8[ . The activation energy is associated with the en- 
ergy required to create a double kink in the vortex which 
is the free energy of two vortex segments along the ab 
plane with the length of a = (Qq/B) 1 / 2 , the average dis- 
tance between the vortices. Using anisotropic Ginzburg- 
Landau theory, the barrier for this plastic movement can 
be estimated as 



U p i — 2E v qq 



d> 2 



87r 2 7 A 2 



B 



1/2 



(2) 



where E v is the vortex energy per unit length along the ab 
planes and A 2 = A a6 A c . With A 2 fc = A 2 b (0)/(1 - 1) (where 
t = T/T c ), an activation energy Uq = U p i oc (1 — t)j\[B 
is predicted which is in qualitative agreement with our 
experimental data. To make an order of magnitude esti- 
mation of Upi, we have used the above value for 7(=13.6) 
and A q6 = 2000 A, thus obtaining U » 1500 K at 1 T, 
which agrees well with the experimental data. 

As the temperature increases, due to the enhanced 
thermal fluctuations, an entanglement-disentanglement 
or 3D-2D decoupling transition should happen. Consid- 
ering the high temperature, we think a 3D-2D decou- 
pling is more likely. Daemen et al. |2Cj have calculated 
self-consistently the decoupling line for a Josephson cou- 
pled superconducting system. When the renormalization 
of the Josephson coupling by thermal fluctuations and 
static disorder are taken into account, the decoupling 
field is given by 



B cr {T) = 



167r 3 /c B TseA 2 h (T)7 2 



(3) 



for moderate anisotropy when £ a {, <C 7s <C A a b, where s 
is the distance between the superconducting Cu-0 planes 
(s - 10 A for Y124) @. Since A 2 b = A 2 b (0)/(1 - t), we 
have B cr (T) oc (1 - T/T c )/T. We find that the B cr (T) 
line we obtained above can be nicely fitted by Eq. (|^). 
The fitted result is B cr = 1406.5(1 - T/T c )/T as shown 
in Fig. 2 by the dashed line. 

As discussed in Ref. 27, the dissipation in the 2D liq- 
uid state is controlled by the plastic motion of pancake 
vortices which can be visualized as the generation of dis- 
location pairs. The typical energy for creating such a 
pair is 



U K = 



le^A^r) 



HBo/B) 



(4) 



with B ~ $ /£ 2 (T). We note that Eq. (§ gives a (1- 
t)lnB dependence of U e as we have just observed. Insert- 
ing the values for s (~10 A ) and A a &(~2000 A ), we find 
U e = 146.5(Zni?o — InB), in excellent agreement with our 
experimental data. The large prefactor po derived from 
TAFF analysis for region h suggests that the liquid in 
region h is very viscous. Since a 2D pancake vortex liq- 
uid is just an extreme case for an entangled liquid, i.e., 
the entanglement length equals to s. The entanglement 
increases the viscosity of vortex liquid and results in a 
larger activation energy, and thus a large prefactor. 

Therefore, our results suggest feasibility of the follow- 
ing scenario: below B m (T), the vortices form a regular 
Abrikosov lattice. As the temperature increases, the lat- 
tice melts into a disentangled liquid, loses its transverse 
correlation while keeping the longitudinal one. The dis- 
sipation is governed by the generation of double kinks. 
Upon further warming, a 3D-2D decoupling transition 
sets in and the vortices lose their longitudinal coherence. 
The 2D vortices form very viscous liquid and the activa- 
tion energy increases which follows a (1 — t)lnB relation- 
ship. 

In conclusion, we have shown that vortex lattice melt- 
ing can also been observed in Y124 single crystals. The 
melting line can be well described by using anisotropic 
GL theory with a Lindermann constant cl = 0.14. Two 
different liquid phases were observed which we think cor- 
respond to the disentangled liquid and the 2D pancake 
vortex liquid, respectively. 
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